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A summary is presented of a scaling study of the finite-temperature chiral phase transition of two-flavor QCD 
with the Kogut-Susskind quark action based on simulations on I 3 x 4 (L—8, 12 and 16) lattices at the quark 
mass of m q = 0.075, 0.0375, 0.02 and 0.01. We find a phase transition to be absent for m q > 0.02, and also quite 
likely at m q — 0.01. The quark mass dependence of susceptibilities is consistent with a second-order transition at 
m q — 0. The exponents, however, deviate from the 0(2) and 0(4) values theoretically expected. 



1. Introduction 

The order of the two-flavor chiral phase transi- 
tion is a basic question in finite-temperature lat- 
tice QCD. Earlier finite-size studies disfavored a 
first-order transition, indicating the chiral transi- 
tion being second-order at a zero quark mass|l|,||. 
This was also corroborated by a scaling analysis 
of Ref. H. To advance the scaling argument, we 
have pushed forward simulations toward larger 
spatial lattice sizes and smaller quark masses, 
where we found Q exponents differing from the 
earlier results We have since completed our 
runsQl, and here present a summary of results. 
For recent similar attempts we refer to Refs. 

2. Simulation 

Our study is made with the plaquette gauge 
action and the Kogut-Susskind quark action on 
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lattices of a size L 3 x 4 with L = 8, 12, 16 at the 
quark masses m q = 0.075,0.0375,0.02,0.01. For 
each set (L, m q ), 10000 trajectories of unit length 
are generated by the hybrid R algorithm for a sin- 
gle value of (3 close to the transition. The stan- 
dard reweighting technique is employed to calcu- 
late the (3 dependence of observables around the 
simulation point. The following susceptibilities 
are calculated at each trajectory: 
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where V — L 3 ■ 4 denotes the lattice volume, Do 
the temporal component of the Kogut-Susskind 
operator, i, j — a, r, and P a , T the spatial and tern- 
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Figure 1. Peak height of x-m and xn as a function 
of spatial volume L 3 . Two runs are made at L = 
12 and m q = 0.01. For xn data plotted with open 
symbols are results from previous studies pip). 



poral plaquette. Averages are taken over the last 
8000 trajectories, with errors estimated through 
jackknife analyses with the bin size of 800 trajec- 
tories. 

3. Finite-size scaling analysis 

In Fig. [j] we plot the results for the spatial vol- 
ume dependence of the peak height of the suscep- 
tibilities Xm and xn- For the latter, results from 
previous studies are also shown for comparison 
(open symbols). For m q > 0.02 the two suscep- 
tibilities do not increase beyond the size L = 12, 
supporting the previous conclusion of an absence 
of phase transition down to m q = 0.025 

For m q — 0.01 an increase of the peak height 
continues up to L = 16 in a linear manner. This 
result is quite different from that of an earlier 
workQ which showed a flattening behavior for 
L = 16 (see open circle for xE ax )- We may as- 
cribe this underestimate to a smaller statistics 
(2500 trajectories^ as compared to 10000 em- 
ployed here). 
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Figure 2. Time history of ipip for the runs with 
m g =0.01. Two runs made for L = 12 are both 
shown. 



A linear increase of susceptibilities is a behavior 
expected for a first-order transition. However, ex- 
amining the time histories of chiral order param- 
eter for m q = 0.01 shown in Fig. ^, we observe 
that metastability signals become weaker toward 
larger spatial sizes: a flip-flop behavior is most 
apparent for L = 8, while irregular fluctuations 
are more dominant for larger L. Correspondingly, 
the histogram of ipip clearly shows a double-peak 
distribution for L = 8, which is less evident for 
L = 12 and hardly visible for L = 16. Further- 
more, if we normalizes L by the pion correlation 
length = l/m w at zero temperature, we find 
the increase of susceptibilities at m q = 0.01 for 
L = 12 — 16 being similar to that at m q — 0.02 
for L = 8 - 12. 

We conclude that the increase of susceptibili- 
ties for m q = 0.01 is probably due to insufficient 
spatial volume, and that a phase transition is also 
likely to be absent for m q = 0.01. 



3 



Table 1 

Critical exponents for each spatial size L as com- 
pared to 0(2), 0(4) and mean-field (MF) values. 
Critical coupling obtained from Xm is employed 
to extract z„. 
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4. Second-order scaling analysis 

For a second-order transition, the critical cou- 
pling and susceptibilities are expected to exhibit 
a scaling behavior toward m q — > given by 
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where the index a labels various susceptibilities 
defined in (1-6). Our results for exponents ob- 
tained by a single-power fit are summarized in 
Table fjj For L = 8 our values are consistent 
with those of Ref. || carried out for L = 8 and 
0.02 < m q < 0.075. 

Th exponent z g determined from the critical 
coupling is roughly consistent with the 0(2) or 
0(4) value theoretically expected. All the other 
exponents exhibit a systematic increase with L 
and deviate significantly from the predictions, 
particularly for zt and z e . 

On the other hand, the hyperscaling relations 
Zg + z m = Zt + 1 and 2zt — z m — z e , which follow 
from the fact that all the exponents are deter- 
mined by the basic thermal and magnetic expo- 
nents, are well satisfied for each spatial size L. 



We have also calculated the scaling function 
F m (x) = m z q m Xm(g 2 ,m q ) with x = {g~ 2 (m q ) - 
g~ 2 (0))-m q z " . We find reasonable scaling if mea- 
sured values of exponents are employed, while re- 
sults are much worse with the use of the 0(4) 
exponents. 

We may summarize that our susceptibility data 
are consistent with a second-order transition, but 
that the exponents apparently take values differ- 
ent from those theoretically expected, at least in 
the range of quark mass m q > 0.01. 

5. Concluding remarks 

This investigation has raised several issues 
which were not apparent in the previous stud- 
ies Jl] [|. While we feel a first-order transition be- 
ing unlikely, a finite-size analysis is needed at a 
larger spatial volume to firmly establish the ab- 
sence of a first-order transition at m q = 0.01. 
The discrepancy of measured values of exponents 
from the theoretical expectations, especially from 
those of 0(2), has to be clarified to confirm the 
second-order nature of the transition. 
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